Zener tunneling in two-dimensional photonic lattices 
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We discuss the interband light tunneling in a two-dimensional periodic photonic structure, as was 
studied recently in experiments for optically- induced photonic lattices [H. Trompeter et al., Phys. 
Rev. Lett. 96, 053903 (2006)]. We identify the Zener tunneling regime at the crossing of two Bloch 
bands, which occurs in a generic case of the Bragg reflection when the Bloch index crosses the edge of 
the irreducible Brillouin zone. Similarly, the higher-order Zener tunneling involves four Bloch bands 
when the Bloch index passes through a high-symmetry point on the edge of the Brillouin zone. We 
derive simple analytical models that describe the tunneling effect, and calculate the corresponding 
tunneling probabilities. 

PACS numbers: 42.25.Bs, 42.82.Et, 42.65.Wi 



I. INTRODUCTION 

Zener tunneling in tilted periodic potentials is an in- 
triguing physical phenomenon which occurs when the en- 
ergy difference imposed on the period by a linear poten- 
tial becomes of the order of the energy gap between two 
nearest Bloch bands A siniple model studied inde- 
pendently by Zener 0, Landau Q and Majorana cap- 
tures the essence of this phenomenon. The well-known 
physical examples include electrical breakdown in Zener 
diodes 0, electrical conduction in nanotubes and su- 
perlattices 'sj, pair tunneling in Josephson junctions 
tunnelling of the Bose-Einstein condensate in optical lat- 
tices Hi Q and an optical analog of tunneling in waveg- 
uide arrays and photonic crystals 0, 0| ■ 

Recent observation of Zener tunneling in two- 
dimensional photonic lattices calls for a nontriv- 
ial generalization of the Landau-Zener-Majorana system, 
since the latter describes only the tunneling at an avoided 
crossing of two Bloch bands. In two dimensions the sit- 
uation of more than two energy levels with comparable 
energy gaps between them is generally unavoidable. This 
is due to the fact that in higher dimensions there exist 
more than one distinct (i.e., not equivalent) Bragg re- 
flection planes. Thus, at least in some cases a n-level 
system with n > 2 must be invoked to describe Zener 
tunneling in a two-dimensional lattice. In this paper we 
analytically study all possible scenarios of Zener tunnel- 
ing in two dimensional square lattices. To this end we 
invoke the so-called shallow lattice approximation (see 
the definition below). Our results can be applied to the 
tunneling in two-dimensional photonic lattices observed 
in Ref. [ij. 

Bloch oscillations and Zener tunneling in two- 
dimensional lattices have been studied numerically in 
Refs. However, a complete analysis of all pos- 

sible cases was not presented. In addition, nonlinear ef- 
fects of Zener tunneling of Bose-Einstein condensate in 
two-dimensional optical lattices has been considered re- 
cently |l4j where it was shown that the modulational 



instability of Bloch waves results in an asymmetry of the 
resonant upper-to-lower band vs. lower-to-upper band 
tunneling (see also Refs. [Tsl [l6l IT^ for other nonlinear 
effects of tunneling in one-dimensional lattices). In our 
case, nonlinearity is negligible and we restrict ourselves 
to the linear Zener tunneling. 

The paper is organized as follows. In section ^ we 
consider the general model for the beam propagation 
in a photonic lattice and discuss the shallow lattice ap- 
proximation. Section IIIII contains the derivation of the 
Landau-Zener-Majorana type models for Zener tunnel- 
ing in a shallow lattice. Conclusions and perspectives 
are summarized in section Hvl 



II. 



MODEL AND SHALLOW-LATTICE 
APPROXIMATION 



We describe the propagation of an optical beam by the 
paraxial equation for the normalized electric field enve- 
lope E 



.dE 1 fd'^E d'^E 



xAn{^)E = 0, (1) 



where ^ — (^1,^2) are dimensionless transverse coordi- 
nates: ^1 = (2Tr/d)x, ^2 = (27r/d)?/; ( = 2ttX/ {nQcP')z 
is the dimensionless propagation distance, with A being 
the wavelength in vacuum, no the bulk refractive index, 
and d the lattice period; x — jnQ{d/ A)^ with 7 being the 
nonlinear coefficient. The total induced refractive index 
pattern An(^) is given as follows 



An 



where the potential 



'cos' 



|]cos^ 



(2) 



(3) 



2 



is due to an optical lattice and the additional potential 
ImiO = y [1 + tanh(Kefl,^)] , 

where k — d/{2TTri), describes the the refractive index 
ramp in the direction e/j (the notations correspond to 
those of Ref. 

A shallow lattice approximation corresponds to the 
condition kAti ^ 1 [see Eq. (^]. For instance, in the ex- 
periment reported in Ref. [ILj the induced optical lattice 
had X = 21.59. Thus to have a weak lattice case within 
the same experimental parameters one should use a weak 
cw laser beam satisfying the condition <C 1/x ^ 0.05. 

Assuming that a weak ramp field (i? ^ 1) changes 
slowly over the lattice period, i.e. 2'7tkB = Bd/rj <C A"^, 
and therefore the Bloch band structure is preserved 0, 
IT^ l20l | , we obtain in this case 

HAn=^ Hlg{^) + ^KB + a^, (4) 

where a ~ kBubbJ^, — jngBd^ /{47rX^ri)eii. Fur- 
ther discussion of validity of the Landau-Zener-Majorana 
models derived below is given in the next section. 
Instead of xlg{$,), we will use the combined potential 

V = Vb[cos(a) +cos(C2) +ecos(a)cos(e2)] (5) 

where Vb = >cA^ /A, and £ is a parameter which enables us 
to establish the effect of the non-separability of the pho- 
tonic lattice Ig{$,) on the tunneling probability. We also 
drop the insignificant constant terms kB/2 and >fA^/4. 
The weak ramp condition now reads \ol\ <C Vq. 
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III. SIMPLE MODELS FOR ZENER 
TUNNELLING 

An optical beam propagating in the lattice experiences 
a strong reflection at the resonant Bragg planes (lines in 
two dimensions), i.e. the Bragg planes which are de- 
fined by the reciprocal lattice vectors for which the lat- 
tice potential has non-zero Fourier components (see, e.g., 
Refs. I^HSl)- Let us first briefly recall the basic theory. 
In the Fourier space the equation for the periodic Bloch 
wave </5q(^), i.e. E — exp{i/3C + i(l$,}'-Pq{^)i reads 

[{q-Qf-p\C^-Q + Y.^Q'-QC^-Q'^Q (6) 
Q' 

where Q denotes a vector of the reciprocal lattice and 
Cq-Q and Vq'_q are the Fourier components of the 
Bloch wave fq{$,) and the lattice y(^), respectively. 
The probe beam with the Bloch index q — {k,£) is 
effectively reflected by the lattice in the first order in 
Vq when the end-point of the Bloch index lies on the 
resonant Bragg plane defined by Q' — Q, i.e. when 
{q - Q? = {q- Q'? = P and Vq-_q ^ 0. 

For the lattice © there are two distinct resonant 
planes, denoted by Bi^q and Bi i in Fig. [21 The cor- 
responding vectors of the reciprocal lattice are Q — 
(2(7b,0) and Q = {2qB,2qB), respectively (in our case, 
qb ~ 1/2). Below we consider the corresponding Bragg 
resonances. 
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FIG. 2: Schematic representation of the extended Brillouin 
zone and the two cases of the Bragg resonance for the two- 
dimensional optical lattice described in the text. Solid lines 
indicate the Bloch index and dashed lines are the resonant 
index given hy q — Q with Q corresponding to the considered 
Bragg plane. 



FIG. 1: Bloch band structure (the first 7 bands) correspond- 
ing to the lattice of Eq. 0. Here Vo = 0.05 and e = 1. 



The lattice (PJ, also called the "quantum antidot" in 
the field of Bose-Einstein condensates j^, corresponds 
to e = 1. Its band structure for Vq = 0.05 is shown in 
figure 



A. Generic Bi o resonance 

The reflection point on the Bragg plane Bi q is q — 
[ibAo) with some £o 7^ ±95 + 0{Vq) (i.e., outside a 
neighborhood of the M-point of radius oc Vb). For a 
broad beam one can proceed in the way similar to the 
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approach of Ref. [23 , that studied accelerating electrons 
in the solid-state theory, and assume that the beam is de- 
scribed by a Bloch wave with the propagation-dependent 
Bloch index q = q{(). By keeping only the resonant 
terms we arrive at the following approximation for the 
Bloch wave describing a broad beam near the point of 
Bragg reflection (or, equivalently, resonance) 

^ = (Ci(C)e^*=(^^«^ + C2(C)e^['('^^-'«"I«^) e*'<^)«^ (7) 

Substituting this expression into Eq. 0J, with kAti = 
V{^) + and projecting on the resonant terms yields 
dk/dC, — — «!, di/dC, = —a2, which are necessary to can- 
cel the linear term in ^ (see, e.g., Ref. 0), and then 
obtain a system of coupled equations for the incident Ci 
and Bragg reflected C2 amplitudes (see also Ref. 0]), 



.dCi 



dC2 



yC2, 



(8) 



[k - 2qB?)C2 



It is 

equal 

aiC- 

exp{- 



convenient to set the propagation variable C 
zero at the resonance point, so that k = qs — 
Redefining the amplitudes as follows (Ci,C2) = 
i{f dC^2 + ql( + a? CV3)/2}(ci, C2) we obtain the 




(9) 



-C2 



Vo 



-Cl. 



Landau-Zener-Majorana system 
. dci 

.dC2 

System ^ is Hamiltonian, its adiabatic energy lev- 
els or, equivalently, the two Bloch bands of (3, /3i,2 = 
Ty/l"^ + (Vb/2)^, experience an avoided crossing at g = 
0, where q = ai(^/2 has the meaning of the running band 
parameter. 

The probability of tunneling P, defined through the 
final amplitude P = |Cj(oo)p for the initial condition 
|Cj(— oo)| = 1, is given by the well-known Landau-Zener 
formula P = exp{— 7rV^^/(2|ai|)} for both the lower- 
to-upper and upper-to-lower band tunneling processes. 
Here ai is the component of a. in the direction perpen- 
dicular to the Bragg plane Bi q, i.e. to the XAf -border of 
the irreducible Brillouin zone. Thus, this case of the in- 
terband tunneling is quasi one-dimensional even it occurs 
in a two-dimensional lattice (e.g., see Ref. (15i..l6..17j1. 



B. Generic resonance 

In this case, the reflection point on the Bragg plane 
Bis is q = (koJo) with ko + £q = 2qB and Iq y^ qs + 
0{Vo). The resonant terms of the Bloch wave read 

■0 = (7^(^)e*'=(f)«i+«^(C)«2_|_(72(C)e*['=(^)^2gB]«i+j[«(C)-2gB]6 

(10)' 



Setting k = ko — ai( and £ = £0 — a2C, redefining the 
amplitudes as (Ci,C2) = exp{—i{[kQ+£Q + {a2 — ai){ko — 
4)]C+[ai+a2]C^/3)/2}(ci, C2) and following all the steps 
described above results in the following Landau-Zener- 
Majorana system 



dci 



-Cl 



-C2, 



(11) 



.dc2 r^iC 



-C2 



-Cl, 



where fli = ai + 02- 

The corresponding Bloch bands are /3i_2 = 
tV?^ + (£Vo/4)2, with q = riiC/2. The interband tun- 
neling probability when crossing the -plane is thus 
given by the formula P — exp{~Tre^VQ / {8\ai + 0:2!)}. 

The Bragg plane -Bi^i is transformed to the TM- 
border of the irreducible Brillouin zone by the follow- 
ing transformation: k = 2qB — k' , £ = £' . Note that 
the transformation changes the coordinates as follows 
^1 = —^'i, ^2 = C2J thus ai — — a'l and 012 = al,. 
Therefore, the probability of tunneling when crossing 
the FM-border of the irreducible Brillouin zone reads 
P = exp{-7re2ynV(8|a2 - ail)}. 

The quantity \a2 — ai\/2 is the component of o: perpen- 
dicular to the FAf-line (exactly \a2 — ai \/V2) multiplied 
by a half length of the respective reciprocal lattice vec- 
tor, \/2qB- In this respect, this case of tunneling is again 
quasi one-dimensional. 



C. Four- fold resonance 

In this more general case, the reflection point is the M- 
point. There exist exactly four Bloch indices which are 
in resonance, they correspond to the four high-symmetry 
points of the first Brillouin zone equivalent to the M- 
point. This case is specific to the two-dimensional tunnel- 
ing, since it leads to the four-level system. The resonant 
part of the Bloch wave reads 

-0 = (:^^(^)e^'=(C)Ci+««(C)?2^,-;2(C)e't''''^^"^''^l^i+'[^('^)"^''^I^= 

_(,C'g((^)g'[fc(C)-29B]Cl-Hi^(C)6 Q(-^)g«fc(C)?l+«[^(C)-29B]C2 

(12) 

Setting k ~ qB — aiC and £ = qB — a2C; defining Cj — 
e«(9BC+[ai+Q2lC /^)Cj and following all the above steps we 
derive the following system 

.dci SliC eVo Va. . 

*rfC " ^Ci + — C2 y(c3 -I-C4), (13) 

.dc2 fliC eVo Vo, . , 

" — ^2 + — Cl + y(c3 -I-C4), (14) 

.dcz n2C sVo Vo. , 

*rfC = - — + — C4 + y(ci+C2), (15) 

.dc4 n2C eVo Vo, . , 

" C3 + y(Cl-t-C2), (16) 
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with f2i = ai + ai and ^i = 0,2 — ol\. The invariance 
of the system with respect to rotations by 7r/2, which 
is the invariance of the lattice itself, is evident from the 
corresponding transformation: a\ — > ai and ai — *■ —a\. 

System (fT^ - (fTB|) can be viewed as a coupling (by the 
terms with Vb /2) of two cases of the Bragg reflection: one 
is for crossing of the _Bi i-plane [Eqs. (|13|l - l|14|) ] and the 
other one is for crossing of the _B_i.i-plane [Eqs. H15|l - 
(|16|l : the reciprocal lattice vector is Q = (— 2(7s, 2(7^))]. 
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FIG. 3: Schematic representation of tunneling at the Af-point 
of the Brillouin zone. The four beams follow the direction of 
the linear ramp (indicated by arrows; here along the FAf-line) . 



A schematic correspondence between the amplitudes 



Cfc 



P of the Bloch wave and the tunneled and reflected 



beams is depicted in Fig. |3| 

The generalized (i.e. multilevel) Landau-Zener- 
Majorana system defined by Eqs. H13|) - (|16|l leads to a 
fourth-order polynomial which does not allow a simple 
analytical solution for the Bloch bands. However, it is 
easy to establish that the four Bloch bands assume the 
following values at the crossing point (C — 0) 



1; 



e 

1'^ 



(17) 



where we have used that \qM\ = l/v^- Thus, recalling 
that in our case e = 1 , we see that the shallow lattice cap- 
tures at least the qualitative behavior of the Bloch bands 
at the Af-point of the experimental two-dimensional op- 
tical lattice of Ref . , where the two inner Bloch bands 
are also crossing at the M-point (see also figure^. 

Though there exists no general formula describing 
the tunneling probabilities between various levels in 
a multilevel Landau-Zener-Majorana system (see, e.g., 
Ref. [23. I2H . l2^ for further discussion), in the generic 
case, i.e. ai_2 7^ 0, the probability of the transition 
between the lowest and the highest adiabatic levels is 
known 0, |2^ . This corresponds to the tunneling across 
the complete band gap discussed in Ref. Assuming 
that \ai + > \a2 — Qfi|, we obtain the probability 
of tunneling between the first and the fourth bands as 
follows 



P — exp < — 



4 Q!2 -I- ai\ 



1 

R 



1 

\a2\ 



(18) 



In the opposite case, i.e. when |ai -|-a2| < |q;2 — the 
Bloch index passes through an equivalent M-point and 
the probability of tunneling is given by the same formula 
with the change \ai + \ ^ \ct2 — cti\. 



1 

0.8 



CM 



0.6 



^ 0.4 



0.2 







-400 -200 







200 400 



FIG. 4: Tunneling at the M-point of the Brillouin zone. Here 
the parameters are: a = 0.05[cos(7r/4), sin(7r/4)], i.e. the 
ramp is in the FA/-direction, Vb = 0.1, and e = 1. Dashed 
line is the analytical prediction given by the formula 1)18^ . 



Two particular cases of tunneling at the Af -point are 
shown in Figs. 0] and El We have simulated numer- 
ically system lfT^ -(|TC )l using the fourth-order Runge- 
Kutta method. Since the system is written for the so- 
called diabatic basis, a large ^-interval was used for re- 
liable results. Asymptotically, i.e. for ( ±00, the 
diabatic basis coincides with the adiabatic one (in our 
case, the latter defines the respective Bloch bands). The 
dashed line gives the analytical prediction of Eq. ifll^ . 
The oscillations of the amplitudes jcfcp in Figs. 01 and 
[31 after the tunneling process reflect a slight difference 
between the respective amplitudes in the diabatic and 
adiabatic bases (in the adiabatic basis there are no oscil- 
lations). 

In Fig. 01 the ramp direction is along the FAf-line and 
the output powers of the side beams (i.e. W oc |C3_4p) 
are the same. In Fig. the ramp direction is close to the 
FAT-bordcr thus creating a substantial asymmetry in the 
output beams 3 and 4, similar to the tunneling observed 
in the experiment pdj . 

We conclude this section with discussion of validity 
of the Landau-Zener-Majorana models considered above. 
In particular, we consider the four-fold _Bi.i-resonance. 
The system (I13f) - (limi approximates the four resonant 
Bloch bands of the lattice at the A/-point (see figure 
The approximate analytical values for the four en- 
ergies at the Af-point, given by equation ((T7jl . can be 
employed to derive a quantitative estimate on validity of 
the Landau-Zener-Majorana model (|13|I - H16|) . The result 
is presented in figure |S1 We have used four values of the 
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FIG. 5: Tunneling at the M-point. Here the parameters 
are: Vo = 0.1, e = 1, and a = 0.05[cos(0.25), sin(0.25)], i.e. 
the ramp is directed close to the FX-line. Dashed line is the 
analytical prediction. 
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FIG. 6: The error of the Landau-Zener-Majorana model for 
the four-fold Bi^i-resonance. The curves give the absolute 
value of the maximal relative difference between the first four 
Bloch bands of the full 2D equation and the corresponding 
energy levels of the four-level Landau-Zener-Majorana model. 



parameter e: e = (solid line; the egg-crate, i.e. sepa- 
rable lattice), £ = 1 (dashed line; the quantum antidot 
in the terminology of Ref. e = 2 (dotted line), and 

£ — —1 (dash-dot line; the quantum dot). 

From figure|n|one can see that for the lattice potentials 
with Vq < 0.05 the relative error of the finite-dimensional 
approximation of Zener tunneling by the Landau-Zener- 
Majorana models is below 5%. For the experimental lat- 
tice of Ref. jOli i-e- the quantum antidot (e = 1), the 
relative error is less then 5% for Vq < 0.15. On the other 
hand, for any particular lattice potential, for the values of 
the lattice strength Vq above some threshold the finite- 
dimensional approximation will break down, as it does 
for the the quantum dot lattice already for Vq = 0.15. 



IV. CONCLUDING REMARKS 

We have derived simple models that allow one 
to describe one-dimensional and more general two- 
dimensional Zener tunneling in two-dimensional periodic 
photonic structures and calculate the corresponding tun- 
neling probabilities. 

We have found that a square two-dimensional photonic 
lattice allows one to observe the quasi one-dimensional 
tunneling, when the Bloch index crosses either the TM- 
border or the XM-border of the irreducible Brillouin 
zone away from the M-point, and the specific two- 
dimensional tunneling, when the Bloch index passes 
through the border of the Brillouin zone at the M-point 
(more precisely, in its neighborhood with radius of the 
order of Vq). We notice that in the first order in Vq the 
Bloch oscillations are not affected by tunneling when the 



Bloch index crosses the FX-border of the irreducible Bril- 
louin zone of the two-dimensional lattice given by Eq. (jSJ 
of section nil 

In our analytical calculations we have employed a 
shallow-lattice approximation which can be related to the 
experimental setup used in Ref. [Tlj only qualitatively, 
and we do observe at least the qualitative correspondence 
with the experiment. There exist two general ways to 
realize a shallow lattice in the experimental setup. The 
first one is to reduce the power of the laser beams creating 
the optical lattice potential Ig{^); this requires the con- 
dition <C l/>f, which we have adopted in this paper. 
The second one is to change the value of the parameter 
X = '^no{d/X)'^ so that >f <C 1. In the latter case, the 
lattice potential will be different from that used in our 
analytic calculations [see Eq. 0], namely it will have 
all higher harmonics in the Fourier expansion. Qualita- 
tively this will result in the following changes. First of 
all, there will be a substantial tunneling when crossing 
the rX-border. Second, the higher harmonics will also 
affect the probabilities of all other tunneling cases. How- 
ever, the latter effect is of a higher order in the potential 
amplitude, since the higher harmonics are not resonant. 
Therefore, the theory developed here can be applied to 
this case as well. The tunneling when crossing the TX- 
border is quasi one-dimensional, and it can be treated 
similar as the quasi one-dimensional tunneling cases dis- 
cussed in this paper. 
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